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ANALYTIC  MODELS  OF  MAGNETIC  FIELD  EVOLUTION  IN 
LASER-PRODUCED  PLASMA  EXPANSIONS 


I.  INTRODUCTION 

A  high  altitude  nuclear  explosion  (HANE)  can  significantly  perturb  the 
natural  ionosphere  and  magnetosphere  by  generating  large-scale  (tens  to 
hundreds  of  kilometers),  long-lived  (several  hours)  ionization  irregularities 
(striations) .  These  irregularities  can  degrade  radar  and  communication 
systems,  e.g.,  through  scintillation  effects.  As  such,  it  is  crucial  to  have 
a  detailed  understanding  of  the  evolution  of  a  HANE  in  order  to  achieve  a 
reliable  communications  system  capability  in  a  nuclear  environment.  Several 
experiments  have  been  conducted^,  both  in  space  (barium  cloud  releases)  and  in 
the  laboratory  (laser-pellet  explosions),  in  order  to  simulate  various  aspects 
of  a  HANE.  Recently,  renewed  interest  in  laboratory  laser-pellet  HANE 
simulation  experiments  has  been  generated  owing  to  detailed  scaling 
studies.*-^ 

An  important  aspect  of  the  early  time  (few  seconds)  evolution  of  a  HANE 
is  the  manner  by  which  the  exploding  debris  plasma  couples  into  the  background 
air  plasma.  The  nature  of  this  early  time  coupling  could  seed  or  influence 
the  evolution  and  structure  of  late  time  ionization  irregularities.  This 
coupling  can  either  by  collislonal  (particle-particle  interactions)  or 
collisionless  (wave-particle  interactions)  depending  on  ambient  densities  and 
temperatures.  Collisionless  coupling  proceeds  via  plasma  microturbulence 
which  in  turn  is  driven  by  various  plasma  instabilities^.  Recently,  a  set  of 
"turn-on"  conditions  for  collisionless  coupling,  in  the  context  of  the  NRL 
laser  plasma  experiment,  has  been  derived^.  A  key  ingredient  in  determining 
whether  or  not  the  aforementioned  plasma  instabilities  will  be  excited  is  the 
structure  and  magnitude  of  the  local  magnetic  field  in  the  debris-air  coupling 
region.  If  the  magnetic  field  compressions  are  too  small  several  of  the 
coupling  Instabilities  will  be  inoperative.  A  detailed  description  of  the 
early  time  magnetic  field  evolution  and  morphology  is  also  important  for 
discussing  related  topics  such  as  electron  heat  transport,  magnetic  field 
driven  interchange  instabilities,  and  ion  leakage  mechanisms.  Analytic  models 
of  the  evolution  of  the  magnetic  fields  can  be  used  to  validate  HANE  numerical 
simulations  and  also  be  compared  with  experimental  results. 
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In  this  report  we  study  analytic  models  of  the  magnetic  field  evolution 
and  compression  in  laser-produced  plasma  expansions.  We  treat  both  one-  and 
two-dimensional  models  both  with  and  without  an  initial  background  magnetic 
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field.  Many  aspects  of  this  problem  have  already  been  investigated  . 
However  these  studies  have  not  been  discussed  or  applied  to  recent  NRL  laser 
pellet  experimental  observations.  The  outline  of  this  report  is  as  follows. 
In  Section  II  we  present  and  discuss  the  general  model  equation  for  the 
evolution  of  the  magnetic  field  in  laser-produced  plasma  expansions.  In 
Section  III  we  study  laser  plasma  expansions  without  an  initial  ambient 
background  magnetic  field  and  the  expected  magnitudes  of  spontaneous  self¬ 
generated  magnetic  fields.  The  calculated  self-generated  magnetic  fields  are 
found  to  be  in  agreement  with  preliminary  NRL  experimental  values.  In  Section 
IV  we  investigate,  using  both  one-  and  two-dimensional  models,  laser  plasma 
expansions  into  an  ambient  background  magnetic  field.  The  predicted  field 
compressions  are  not  inconsistent  with  those  obtained  using  NRL  experimental 
observations.  We  find  that  one-dimensional  models  using  sharp  debris  density 
profiles  give  unrealistically  large  magnetic  field  compressions.  We  show  that 
two  dimensional  models  with  diffuse  profiles  can  explain  several  experimental 
observations.  Finally  in  Section  V  we  summarize  our  results. 

II.  BASIC  EQUATIONS 

2, 

For  high  B  (B  /8rr  «  Nk^T)  plasmas,  the  equation  for  the  evolution  of  the 

B  I?-!  A 

magnetic  field  B_  can  be  written1^-1^ 

3B_  2 

—  -VxVxB-^Vx[ri»(VxB)]+S  (1) 

where 

ckB 

S  -  -  VN  x  VT  (2) 

—  eN  e  e  v  ' 

e 

with  _n  the  resistivity  tensor,  V_  the  fluid  velocity,  c  the  speed  of  light,  k^ 
is  Boltzmann's  constant,  and  Ng  and  Tg  are  the  electron  temperature  and 
density,  respectively.  Radiation  pressure  effects  have  been  neglected  in  Eq. 
(2)  since  they  can  be  shown  to  be  small  for  the  laser  intensities  used  in  the 
current  NRL  experiment.  Equation  (1)  is  simply  Faraday's  law  with  the 
electric  field  determined  from  the  force  balance  equation  for  the  electrons. 


2 


The  generation  of  a  magnetic  field  requires  that  the  last  term  in  Eq.  (1),  the 
source  term  J_,  be  nonzero.  This  requires  that  7Ng  and  7Tg  be  nonparallel. 

By  defining  the  dimensionless  quantities  t  =  VQt/L,  \f  ■  V/Vq,  i  “  x/L, 

Eq.  (1)  can  be  written 


-  *  7  x  if  x  B  -  r'1  f! 

v  — -  m  — 

at 


-1  -2 
>  *  f*  ^  f 


(3) 


where,  for  the  moment,  we  have  neglected  _S_  and  defined  VQ  and  L  as  a 
representative  fluid  velocity  and  a  magnetic  gradient  scale  length, 


Here  Rm  defines  a  magnetic  Reynolds  number  Rm=  4naLVQ/c2 


respectively. 

with  n=I  o  *  and  I  the  unit  tensor.  For 


>  1,  magnetic  field  convection 


(V  x  xB^  term)  dominates  over  diffusion  (v  B  term)  whereas  for  Rm  <  1  the 


opposite  is  true. 


eff 


An  effective,  electron  collision  frequency  v  can  be 


eff 

defined  by  o  =N  e2/m  v  where  N  and  m  are  the  electron  density  and  mass. 


respectively.  For  v 


Iff 


ei 


.  e  -  “e 
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3  x  10  Z  £n  A  N  /T  sec  ,  the  classical 
e  e 


Coulomb  collision  frequency,  with  Z  the  charge  number  and  £n  A  the  Coulomb 
logarithm,  we  find  R  =  104  L(cm)  where  we  have  taken  T  a  100  eV  and 


m  '  '  e 

V  a  4  x  107  cm/sec  as  representative  NRL  laser  experimental  parameters  (B. 


Ripin,  private  communication).  For  L  ~1  cm,  Rm  a  104  dissipative  effects  are 
negligible.  However,  if  the  effective  collision  frequency  in  increased  by 
other  processes,  e.g.,  plasma  microturbulence,  Rjjj  will  decrease  and  resistive 
effects  will  become  more  important.  For  example,  for  plasma  turbulence 
resulting  from  the  magnetized  ion-ion  instability^  an  effective  collision 
frequency  can  be  written 


eff 


'ij 


n  1S  Pi  (  2/3  ,-1/3  ,1/2  1/3  2/3. % 

°*15  “Hi  7"  (“ji  +  2  3  (“ji  “ji 
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ji 


‘ji 


(4) 
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where  i  and  j  refer  to  ion  species  i  and  j,  ■  upi^  +  » 


pe 


pe 


and  electron  plasma  frequencies,  is  the  mass  density  of  species  i,p  the 


total  mass  density  and  a,  .  "  N,  Z  2  m. /N,  Z  2  m,.  For  aluminum  (i)  streaming 

1 J  _3 

through  nitrogen  (j),  i.e,  a  101®  cm  ,  a  1014  cm  , 


Zt  *  10,  Z4  *  3  we  find  Rffl  a  10  L(cm)  -  10  for  L  a  1  cm  making  resistive 


i  ’  j 
effects  more  important 


III.  LASER  PLASMA  EXPANSION  WITHOUT  A  BACKGROUND  MAGNETIC  FIELD 


As  shown  in  the  previous  section,  spontaneous  magnetic  fields  can  be 
generated  by  nonparallel  density  and  temperature  gradients.  Thene  self¬ 
generated  magnetic  fields  can  be  quite  large*^  near  the  focal  spot  region. 
Since  these  spontaneous  fields  will  be  carried  along  with  the  expanding  plasma 
they  could  influence  greatly  the  electron  and  ion  dynamics  in  the  coupling 
shell.  In  addition,  large-self  fields  also  imply  asymmetric  departures  from 
completely  spherical  expansion  and  reduced  coupling  of  the  laser  engergy  into 
the  target. 

According  to  Eq.  (2),  the  magnitude  and  direction  of  the  self  generated 

magnetic  field  is  determined  by  the  geometrical  configuration  of  the  laser- 

plasma.  A  laser  beam,  which  is  cylindrically  symmetric,  will  produce  a  plasma 

which  expands  in  the  direction  of  the  normal  to  the  target  and  is  symmetric 

about  its  expansion  direction.  From  symmetry  considerations  there  can  be  no 

azimuthal  density  or  temperature  gradient.  During  the  laser  heating  of  the 

target,  it  is  reasonable  to  assume  that  the  largest  contribution  to  the  source 

term  in  Eq.  (2)  comes  from  a  temperature  gradient  in  the  radial  direction  and 

a  density  gradient  in  the  direction  of  the  target  normal  due  to  expansion  of 

the  target  plasma.  Due  to  the  finite  radial  extent  of  the  laser  beam  a  radial 

temperature  gradient  will  exist  near  the  edge  of  the  focal  spot.  This 

combination  of  V  T  and  7  N  will  generate  a  magnetic  field  in  the  azimuthal 
6  6  12 

direction  in  the  form  of  a  torus  .  The  self-generated  field  will  be 
convected  radially  by  the  expanding  plasma.  This  scenario  has  been  confirmed 
by  many  previous  laser-pellet  experiments*^. 

In  order  to  approximate  the  self-generated  magnetic  fields  in  the  NRL 
laser  HANE  experiment  we  assume  a  purely  radial  temperature 
gradient  3Te/3r  and  a  density  gradient  3Ng/3z  with  r  denoting  distance 
perpendicular  to  the  normal  to  the  target  plane  and  z  representing  distance  in 
the  axial  direction  perpendicular  to  the  target  plane.  As  a  result  Eq.  (2) 
gives 

S  -  (3B/3t)self  “  (<VeNe)  (3Ne/3z)  (3Var)  <5> 

The  radial  temperature  gradient  is  of  the  order  3T  / 3r  «  T  /r  where  r  is  the 

e  e  o  o 

radial  extension  of  the  laser  heated  plasma  near  the  focal  spot.  We  take  the 


density  gradient  in  the  z-direction  to  be  given  by  the  debris  ion  expansion 
velocity  V  ,  i.e,  3lnN^/3z  *  3lnNg/3z  =  [V0tl+  (AVo/Vq)R]  ^  where  xL  is  the 
duration  of  the  laser  pulse,  AVq  is  the  thermal  debris  velocity  spread,  and  R 
is  the  approximate  position  of  the  debris  density  maximum.  Taking 
(3B/3t)gelf  *  Bself/TL  we  have  from  (5) 

®aelf  ■  9  *  107  ITT  t1  +  <4W  (R/VotL>]'1  G  (6) 

o  o 

where  Te,  VQ,  rQ,  x^  are  expressed  in  eV,  cm/sec,  cm,  and  sec, 

respectively.  For  T  *100  eV,  V  *4  x  107  cm/sec,  r  *lcm,  AV  /V  =  0.2, 
v  3  e  _9  0  o*oo  * 

R  *  0.5  cm,  tl  ■  4  x  10  sec  (B .  Ripin  private  communication),  Eq.  (6) 

gives  B  j*  *  100  G  which  is  in  agreement  with  experimentally  measured  values 

(S.  Kacenjar,  private  communication). 

To  find  the  approximate  time  dependence  of  the  self-generated  magnetic 
1  o 

field,  we  consider  the  fluid  variables  Ne,  Te,  V  in  Eq.  (1)  as  consisting  of 

zeroth  order  contribution  plus  a  first  order  part,  i.e., 

N  ■  N  +  AN  ,  T  *  T  +  AT  ,  and  V  *  V  +  AV  .  The  zeroth  order  parts  N__, 
e  eo  e  e  eo  e  o  r  eo* 

Teo,  and  VQ  describe  a  spherically  symmetric  expansion  with  the 

perturbations  A  N  ,  A  T  ,  and  A V  representing  a  small  departure  from  spherical 

symmetry  giving  rise  to  a  source  term  vanishes  in  the  spherically  symmetric 

case.  In  other  words,  we  can  linearize  Eq.  (1)  and  solve  for  3AB/3t  using  the 

zeroth  order  motion  for  V_  and  in  the  first  two  terms  on  the  right  hand  side 

of  Eq.  (1).  The  perturbations  A  N  and  A  T  are  retained  in  the  source 

e  e 

term  j>. 

The  radius  of  the  expanding  laser  plasma,  rg,  is  found  from^  the 
following 

s  (t)  ■  R»  +  il  \  (t')  dt'  <7> 

with  Rq  the  initial  radius.  The  velocity  of  the  expanding  plasma,  in 
cylindrical  coordinates,  follows  from  the  continuity  equation 

Vr  “7  Vs  <C>  <8a> 

s 
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The  exploding  plasma  is  assumed,  In  zeroth  order,  to  have  a  homogeneous 
density  N  (t)  in  the  course  of  its  assumed  adiabatic  expansion  so  that 


I 

T 


3(y-l) 


(9) 


Furthermore,  we  assume  that  the  expansion  decreases  the  plasma  temperature  so 
that  the  conductivity  scales  as 


a 


£  )3/2 

o 


9(y-l)/2 


(10) 


where  oq  is  the  initial  value  and  a  Coulomb  collisional  conductivity 


15 


has 


been  assumed.  As  a  result,  Eq.  (l)-(2)  for  the  azimuthal  component  A  B  of 
the  self-generated  field,  asuming  cylindrical  symmetry  (—  **  0)  can  be  written 

3Ba  a  a  l  a  3Ba  Bs  a2Ba 

(V  B  )  -  (V  B  )  +  Dfi  r  - - £  +  -r— it— ) 

3t  3z  z  8  3r  r  o'  '■r  3r  3r  rz  3zz  1 


+  S  (r,  z,  t) 


with  S 


ck_  3N  3T  3N  3T 

o  f _ e  _ e _ e  _ e 

eN  lar  3z  3z  3r  ' 

e 


do 


where,  for  clarity,  we  have  dropped  the  A"s  from  BQ,  V,  Ne>and  Tg  and 

D(t)  ■  D  (r  /R  )^Y  with  D  -  cz/4iro  We  with  to  solve  Eq.  (11)  as  an 
o  s  o  o  o 

initial  value  problem  with  B„  (t**0)  «0  and  V_,  V_  given  by  Eq.  (8). 

0  £  Z 

1 2 

For  the  expanding  plasma  cloud,  we  consider  the  following  debris 
density  and  temperature  profiles  which  are  smooth  functions  of  position. 


Ne  "  Neo  “A  N  e  exp  (~r2/r2s(t)) 

T  "  "  (AT  /r  )  /*  dz-  exp  (-z'z/rz  ) 

eeo  es'o  s 

From  (12a)  and  (12b)  we  find 


S  (r,  z,  t)  -  (-2  c  kjj/e)  (A  Ne/Ng)  (r/rs)  (A  Tg/rz) 


(12a) 

(12b) 


x  exp  [-(r  +  z  )/rg] 


(13) 


Tidman  (1975)  has  solved  Eq.  (11)  using  the  profiles  from  (12)  and  finds,  in 
scaled  time  units, 

,  „  AN  AT  V 

B  (t)  -  -  8w  (2  e2)-1/2  Bq  (-^)  (^-)  (/)  t  (1  +  t)  (14) 

e  e  o 

where  B  ■  k„  T  o  /ec  ,V,=D/R  =c2/4ito  R  ,r  =  R  +  V  t  (constant 
oboo  doo  oosoo 

expansion  velocity  Vq),  t  3  Vq  t/RQ  and  e  =  2.71828  .  Using  parameters 

typical  of  the  NRL  laser  experiment  (Tg  =  100  eV,  Ng  =  10^  cm~\  VQ  = 

4  x  102  cm/sec),  B  (t)  is  given  by  Fig.  1  for  several  ilues 

0 

of  A  N  /  N  and  AT  /  T  .  Both  the  magnitudes  (several  hundred  and 

e  e  e  e 

time  dependence  of  B  are  consistent  with  preliminary  measurements  •  zed 

0 

from  the  current  NRL  laser  HANE  experiment. 

IV.  LASER  PLASMA  EXPANSION  WITH  A  BACKGROUND  MAGNETIC  FIELD 

For  realistic  simulations  of  a  HANE,  a  background  magnetic  field  must  be 
introduced  into  laser-pellet  experiments.  With  a  background  magnetic  field, 
magnetic  field  compression  can  now  take  place  in  addition  to  spontaneous 

magnetic  field  creation  as  disucssed  in  Section  III.  Magnetic  field 
compression  may  be  the  first  stage  in  the  process  leading  to  "pickup"  of  the 
background  air  ions.  It  is  important  to  compute  the  spatial  and  temporal 
history  of  the  magnetic  field  compression  in  order  to  determine  where  and  when 
the  peak  compression  is  achieved. 

A.  One-dimensional  models 

We  consider  a  one-dimensional  model  of  laser-induced  plasma  expansions 
into  a  background  magnetic  field.  The  model  consists  of  a  debris  plasma 

a 

streaming  with  velocity  through  a  stationary  background  (air)  plasma. 

Choosing  to  work  in  the  debris  frame  of  reference,  the  ion  component  of  the 
expanding  debris  plasma  is  stationary  with  density  nD  while  the  background 

A 

plasma  is  assumed  uniform  with  density  ng  and  having  flow  velocity  -V^jc  •  The 

basic  configuration  is  shown  in  Fig.  2.  In  the  interaction  region 

-x  <  x  <  0  ,  continuity  and  quasi  neutrality  are  imposed  where  x  =  V  tt 

O  O  U  L 

with  tt  the  duration  of  the  laser  pulse.  These  conditions  determine  the 
u  - 

density  ng  ■  n^  +  tip  and  flow  velocity  -  x  ■  -  (n^/n^  +  n^)  x  of  the 
debris  electrons.  Initially,  a  constant  background  magnetic  field  B  z  is 


taken  to  be  normal  to  the  flow  but  excluded  from  the  Interaction  region  as 
shown  in  Fig.  2.  We  wish  to  determine  the  field  compression  Bz  (x,t)  /  Bq  . 

For  the  case  where  collisions  are  absent,  the  evolution  of  (x,t)  was 

determined  by  Longmire^  using  magnetic  flux  conservation  arguments  and  is 

illustrated  in  Fig. 3.  Here  Bz  jumps  discontinuously  from  BQ  at 

X  =  e  to  (V.  /V)B  >  B  atX=-£.  The  leading  edge  of  the  compression  in 
d  e  o  o 

the  interation  region  is  convected  with  velocity  Vg. 

Including  collisional  effects  Eq.  (1)  gives  for  Bz  (x,t) 


—  (V  B  )  +  -5-  V  (x)  L  (x)  B 
3x  x  z  3x  x  c  3x  z 


(15) 


where  L  (x)  =  v  (x)  c2/io  2  (x)  V  (x)  with  the  effective  collision 
c  pe  x  _  ^ 

frequency  v  (x)  being  defined  from  o  (x)  =  v  (x)  mg  /  N  (x)  c2  .  It  should 
be  noted  that  in  the  derivation  of  Eq.  (1)  and  (2)  the  effects  of  electron 


inertia  have  been  neglected. 


cases 


7,8,10 


Eq.  (15)  has  been  solved  for  several  special 


As  a  simplified  model  we  consider  the  case  where  the  effective  collision 
frequency  is  non  zero  only  in  the  interation  region,  -xq  <  x  <  0  .  This 

effective  collision  frequency  is  defined  as  the  sum  of  t)  •»  classical  Coulomb 

collision  frequency  plus  an  anomalous  part  due  to  plasma  turbulence.  Let 

0  x  >0 

L  (x)  =  {  L  -x  <  x  <  0 

c  1  co  o 

0  x  <  -x 

where  Lco  is  constant.  As  a  result  the  equation  governing  Bz(x,t) 

for  -x  <  x  <  0  can  be  written 
o 


0 


(16) 


In  the  following  dimensionless  variables  x  ■  x  /  L 
and  B  *  Bz/®0  Eq.  (16)  can  be  written  (  -x  <  x  <  0) 


Vllt  1  Lco  • 


where,  for  clarity,  we  have  dropped  the  tildes.  The  initial  and  boundary 
conditions  appropriate  to  Eq.  (17)  are:  B(x  >  0.  t  ■  0)  ■  1,  B(-x  < 


x  <  0,  t  -  0)  -  0,  B(x  -  «,t)  -  1,  B(x“0-e,t)  +  3B/3x  (x-0-e,t)  “  l. 

Eq.  (17)  together  with  these  Initial  and  boundary  conditions  can  be  solved  to 


yield 


<4  -1 

B  (-  x  <  x  <  0,  t)  -  tt—  (  1  +  (2  w)  exp  (-  x/2)  j;(x.t)) 

Z  O  V  v  ' 


where 


5  (x,  t)  -  (1  -  2  3  /  3  x)  a  (x,  t) 


CD  \  -2  1 

a  (x,  t)  ■  J_  d  u  u  (u2  +  -r  )  exp  (  -  (u2  +  -r  )t  )  sin  u  x 


-  tr  (x  -  t)  exp  (-|x  1/2)  +  (it  /2) 


ft  /  -  3/2  -1/2..  f-x2  m-. 

/  d  m  (  t  m  -m  )  exp  (-r—  -  j) 


A  u  t  *■  “  I4m  4 1 

m  x  (x-t)exp(-|x| /2)  +  nerfc(i  x  t  ^2)  Z  x-  -  —n+t ^ - 

n-o  2  (2n-l)!!n! 

i/o  ”  2n,.  . ^ \  n-1  ...  .  k+l/2 

+  x1/2exp(-x2/4t)  Z  Z  (-2)k  1(2k-l)!!  (^-) 

n-l  2  (2n-l)!!n!  k-o  x 


Fig.  4  shows  the  evolution  of  Bz(x,t)  as  given  by  Eq.  (18).  For  times 
t  >  1  collisions  will  smooth  the  discontiuous  jump  In  Bz  given  the  Longmire 
analysis  • 

B.  Two-dimensional  models 

In  the  previous  section,  a  one-dimensional  model  of  an  expanding  laser 
plasma  was  assumed  to  have  sharp  discontinuous  debris  density  profile  In  order 
to  simplify  analytical  calculations.  Unrealistically  large  magnetic  field 
compressions  were  obtained^  by  assuming  an  Infinitely  steep  density  profile. 
In  this  section  we  show  the  effects  of  using  a  diffuse,  smooth  debris  density 
profile  by  making  a  two-dimensional  analysis^  of  a  laser  plasma  expanding 
Into  a  background  magnetic  field. 


consider  a  laser  plasma  expanding  into  a  paraxial  magnetic 


B  z 
o 


B  cos  0  e  -  B  sin  9  eQ  which  is  imbedded  in  a  stationary 
o  r  o  9 


We 

field  B 

background  plasma  and  assume  a  spherical  coordinate  system  (r,  9,  $  )  with  the 
polar  angle  9  measured  from  the  z-axis.  Assuming  spherical  symmetry 
(  3/3  9*3/  3  $  ■  0  )  expressions  for  Br  and  B^  from  Eq.  (1)  can  be  written 

3B_ 


ST  ‘  (  '  ta"  6)_1  Vr  B9 


„  „  3B  2B  2Bft  „  . 

+  D  (r-2  4  -ji  -  — g-*— — — — — •  ) 

^  3r  3r  r2  r2  1 


3B , 


3t 


-  r  _1  r  v  Bo 

3r  r  9 


+  D  ( 


-2  3 


3B 


9 


9 


3r 


3r  r2  sin2  9 


(19) 

(20) 


where  we  have  neglected  the  self-generated  magnetic  fields  in  comparison  with 

Bn  and  noted  the  B  is  negligible  for  very  early  times.  On  time  scales  short 
u  <p 

compared  to  the  resistive  time  scale,  the  diffusive  terms  in  (19)  and  (20)  can 
be  neglected  and  we  are  left  with 


I 


3B 

3t 

3B 

3t 


£ 

9 


(r  tan  9)  ~1  BQ 


V 

r  9 


(21) 

(22) 


We  assume  uniform  expansion  vr”vor^rO*  9uas*-’neutral:I-ty  f°r  the  electrons,  and 
a  gaussian  debris  density  profile  ND(r,t-0)-N0exp  (-r2/r2)  where  r0*voxL 
with  tl  the  duration  of  the  laser  pulse.  Under  these  conditions,  Eqs  (21)  and 
(22)  can  be  solved11  analytically  for  B^,  Br  for  short  times  6t  such  that  the 
field  is  only  slightly  altered: 


„  2  VQ  Be  at  (No/Nb  )  exp  (-  r2  /  rQ2) 

®9  ro  f-1  +  <No  /  6XP  rZ  ^  ro 


N  22 

-  exp  (~r  /ro) 


SB 


r 


(r  tan  9) 


N  V 


t  o  o 
lNB  +N0 


(23) 


(24) 


where  is  the  background  plasma  density.  The  peak  magnetic  compression  is 
achieved  at  r*»rm  where  3  /  3  r  (6  B^  )  ■  0  giving 


(25) 


<  s  B„  1  Bo  ’max  *  <Vo  S  ‘  >  2  ro  >  (  (  rm  1  ro  >  '  2  ) 

1/2 

with  r  =  (  f  (N  ,  fL.  ))  r  ,  f(x,y)  »  ln(x/y)-ln  g(x,y),  g(x,y)= 

(An  (x  /  y)  +  2)  /  (in  (x/y)  -  2)  .  For  example  in  1L  =  1014  cm  and 

-3  B 

N  =  101®  cm  ,  we  find  r  =  1.9  r  Indicating  that  peak  magnetic  field 
o  mo 

compression  occurs  ahead  of  the  expanding  laser-produced  plasma  cloud. 
Furthermore,  the  laser  plasma  cloud  density  at  rm  is  given  by 

N  (r  =  r  )=  0.28  N„  .  A  characteristic  shell  thickness,  at  early  times,  can 
o  m  B  i 

be  found  using  6  B  =  -*■  (6  B  )  ,  giving  t  =  0.42  r  .  These  analytically 

Z  o  UlclX  O 

obtained  features  that  (1)  the  maximum  field  compression  occurs  in  front  of 
the  advancing  shell  and  (2)  the  maximum  field  intensity  in  the  compression 
region  is  proportional  to  the  radial  displacement  of  the  shell  are  consistent 
with  preliminary  NRL  laser  experimental  results  (S.  Kacenjar,  private 

communication) . 

Finally,  by  consideration  of  the  conservation  of  magnetic  flux  in  two- 
dimensions,  the  field  compression  scaling  as  given  by  Eq.  (23)  and  (24)  can  be 
shown  to  be  reasonable.  Consider  an  annulus  of  compressed  magnetic  field  with 
inner  and  outer  radii  of  and  respectively.  Magnetic  flux  conservation 
implies  ®or2  =®C(r2  “rl  '  where  Bc/Bq  is  the  magnitude  of  the  field 

compression.  This  gives  Bq/  Bq  =  r^  /  2  6  *  ®  “  r2  ~  rl  '  For  the  NRL  laser 

experiment  R2  =  3cm,  6  =  1.3cm  ,  giving  Bc/Bq  -  1“2  in  agreement  with 

observations  (S.  Kacenjar,  private  communication).  In  addition,  the  direct 

proportionality  between  B£/Bq  and  r£  is  not  inconsistent  with  recent  NRL 

experimental  findings  (S.  Kacenjar,  private  communication). 

V.  SUMMARY 

We  have  attempted  to  provide,  in  this  preliminary  report,  simple  analytic 
models  for  the  magnetic  field  evolution  in  laser-induced  plasma  expansion  in 
the  context  of  the  NRL  laser  plasma  HANE  simulation  experiment.  Both  one-  and 
two-dimensional  models  have  been  used  for  laser  plasma  expansions  with  and 
without  initial  background  magnetic  fields.  For  the  case  with  no  initial 

background  magnetic  field,  we  find  reasonable  agreement  between  analytic  two- 
dimensional  models  and  preliminary  results  from  the  NRL  laser  plasma 
experiment.  It  is  shown  that  self-generated  spontaneous  magnetic  fields  are 
small  in  coraparision  to  proposed  ambient  background  fields.  For  the  case  of 


,3 


Fig.  1  Self-generated  magnetic  field  Ba  produced  by  an  asymmetlc  laser- 

0 

plasma  expansion  vs.  scaled  time  t  -  Vot/RQ.  Curve  a  corresponds 
to  AN/N  ■  AT/T  -  0.1  while  curve  b  corresponds  to  AN/N  -  AT/T  *  0.3. 


Fig.  2  Plot  of  (a)  density  of  debris  Ions  N0,  background  Ions  Ng,  and 
electrons  Ng  ■  ND  +  Ng  as  a  function  of  x  and  (b)  Initial  (t  -  0) 
profile  of  background  magnetic  field.  Here  xQ  ■  Vdtl  with  VD  the 
average  debris  velocity  and  tt  the  duration  of  the  laser  pulse. 


Plot  of  Bz(x,t)  with  constant  resistivity  for  -  xQ  <  x  <  0.  T 
quantities  x,  t,  Bz  are  scaled  by  LCQ,  Lco/Ve,  and  (VD/Ve)B 
respectively. 
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WASHINGTON,  D.C.  20330 
01CY  ATTN  AFRDQ 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01731 
01CY  ATTN  J.  DEAS 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION/YSEA 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01732 
OICY  ATTN  YSEA 


HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION/DC 
DEPARTMENT  OF  THE  AIR  FORCE 
ilANSCOM  AFB,  MA  01731 

01CY  ATTN  DCKC  MAJ  J.C.  CLARK 

COMMANDER 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
WRIGHT-PATTERSON  AFB,  OH  45433 
01CY  ATTN  NICD  LIBRARY 

01CY  ATTN  ETDP  B.  BALLARD 

O  MMANDER 

RO.-Ll-,  AIR  DEVELOPMENT  CENTER,  AFSC 
GRir’FISS  AFB,  NY  13441 

0  ICY  ATTN  DOC  LIBRARY/TSLD 
0  ICY  ATTN  OCSE  V.  COYNE 

SAMSO/SZ 

POST  OFFICE  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 
(SPACE  DEFENSE  SYSTEMS) 

01CY  ATTN  SZJ 

STRATEGIC  AIR  COMMAND/XPFS 
OFFUTT  AFB,  NB  68113 

01CY  ATTN  ADWATE  MAJ  BRUCE  BAUER 
0 ICY  ATTN  NRT 

0  ICY  ATTN  DOK  CHIEF  SCIENTIST 

SAMSO/SK 
P.O.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 

01CY  ATTN  SKA  (SPACE  COMM  SYSTEMS) 
M.  CLAVIN 

SAMSO/MN 

NORTON  AFB,  CA  92409 
(MINUTEMAN) 

01CY  ATTN  MNNL 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB,  MA  01731 

01CY  ATTN  EEP  A.  LORENTZEN 

DEPARTMENT  OF  ENERGY 
LIBRARY  ROOM  G-042 
WASHINGTON,  D.C.  20545 

01CY  ATTN  DOC  CON  FOR  A.  LABOWITZ 


DEPARTMENT  OF  ENERGY 
ALBUQUERQUE  OPERATIONS  OFFICE 
P.O.  BOX  5400 
ALBUQUERQUE,  NM  87115 

01CY  ATTN  DOC  CON  FOR  D.  SHERWOOD 

EGS.G,  INC. 

LOS  ALAMOS  DIVISION 

P.O.  BOX  809 

LOS  ALAhoS,  NM  85544 

01CY  ATTN  DOC  CON  FOR  J.  BREEDLOVE 

UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  808 
LIVERMORE,  CA  94550 

01CY  ATTN  DOC  CON  FOR  TECH  INFO  DEPT 
01CY  ATTN  DOC  CON  FOR  L-389  R.  OTT 
01CY  ATTN  DOC  CON  FOR  L-31  R.  HAGER 
01CY  ATTN  DOC  CON  FOR  L-46  F.  SEWARD 

LOS  ALAMOS  NATIONAL  LABORATORY 

P.O.  BOX  1663 

LOS  ALAMOS,  NM  87545 


01CY 

ATTN 

DOC 

CON 

FOR 

J. 

WOLCOTT 

01CY 

ATTN 

DOC 

CON 

FOR 

R. 

F.  TASCHEK 

01CY 

ATTN 

DOC 

CON 

FOR 

E. 

JONES 

01CY 

ATTN 

DOC 

CON 

FOR 

J. 

MALIK 

01CY 

ATTN 

DOC 

CON 

FOR 

R. 

JEFFRIES 

01CY 

ATTN 

DOC 

CON 

FOR 

J. 

ZINN 

01CY 

ATTN 

DOC 

CON 

FOR 

P. 

KEATON 

01CY 

ATTN 

DOC 

CON 

FOR 

D. 

WESTERVELT 

01CY 

ATTN 

D.  SAPPENFIELD 

SANDIA  LABORATORIES 
P.O.  BOX  5800 


ALBUQUERQUE, 

NM 

87115 

01CY 

ATTN 

DOC 

CON 

FOR 

W.  BROWN 

01CY 

ATTN 

DOC 

CON 

FOR 

A.  THORNBROUGH 

01CY 

ATTN 

DOC 

CON 

FOR 

T.  WRIGHT 

01CY 

ATTN 

DOC 

CON 

FOR 

D.  DAHLGREN 

01CY 

ATTN 

DOC 

CON 

FOR 

3141 

01CY 

ATTN 

DOC 

CON 

FOR 

SPACE  PROJECT  1 

SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE,  CA  94550 

01CY  ATTN  DOC  CON  FOR  B.  MURPHF.Y 
01CY  ATTN  DOC  CON  FOR  T.  COOK 

OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  D.C.  20545 

01CY  ATTN  DOC  CON  DR.  YO  SONG 
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OTHER  GOVERNMENT 

DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON,  D.C.  20234 

01CY  (ALL  CORRES:  ATTN  SEC  OFFICER  FOR) 

INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMUNICATIONS  &  INFO  ADMIN 
BOULDER,  CO  80303 

01CY  ATTN  A.  JEAN  (UNCLASS  ONLY) 

01CY  ATTN  W.  UTLAUT 
01CY  ATTN  D.  CROMBIE 
01CY  ATTN  L.  BERRY 

NATIONAL  OCEANIC  &  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 
01CY  ATTN  R.  GRUBB 
01CY  ATTN  AERONOMY  LAB  G.  REID 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
01CY  ATTN  I.  GARFUNKEL 
01 CY  ATTN  T.  SALMI 
01CY  ATTN  V.  JOSEPHSON 
01CY  ATTN  S.  BOWER 
01CY  ATTN  D.  OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

01CY  ATTN  RADIO  SCIENCES 

AUSTIN  RESEARCH  ASSOC.,  INC. 

1901  RUTLAND  DRIVE 
AUSTIN,  TX  78758 
01CY  ATTN  L.  SLOAN 
01CY  ATTN  R.  THOMPSON 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  BOX  983 
BERKELEY,  CA  94701 
01CY  ATTN  J.  WORKMAN 
01CY  ATTN  C.  PRETTIE 
01CY  ATTN  S.  BRECHT 


BOEING  COMPANY,  THE 
P.O.  BOX  3707 
SEATTLE,  WA  98124 
01CY  ATTN  G.  KEISTER 
01CY  ATTN  D.  MURRAY 
01CY  ATTN  G.  HALL 
01CY  ATTN  J.  KENNEY 

CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 
01CY  ATTN  D.B.  COX 
01CY  ATTN  J.P.  GILMORE 

COMSAT  LABORATORIES 
LINTHICUM  ROAD 
CLARKSBURG,  MD  20734 
01CY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

01CY  ATTN  D.T.  FARLEY,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 
01CY  ATTN  H.  LOGSTON 
01CY  ATTN  SECURITY  (PAUL  PHILLIPS) 

EOS  TECHNOLOGIES,  INC. 

606  Wilshire  Blvd. 

Santa  Monica,  Calif  90401 
01CY  ATTN  C.B.  GABBARD 

ESL,  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 
01CY  ATTN  J.  ROBERTS 
01CY  ATTN  JAMES  MARSHALL 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 

01CY  ATTN  M.H.  BORTNER  SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.O.  BOX  1122 
SYRACUSE,  NY  13201 
01CY  ATTN  F.  REIBERT 
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GENERAL  ELECTRIC  TECH  SERVICES  CO.,  INC. 
HMES 

COURT  STREET 
SYRACUSE,  NY  13201 

01CY  ATTN  G.  MILLMAN 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 
01CY  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
01CY  ATTN  TECHNICAL  LIBRARY 
01CY  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

GTE  SYLVANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

01CY  ATTN  DICK  STEINHOF 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

01CY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 

OICY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA,  VA  22311 
OICY  ATTN  J.M.  AEIN 
OICY  ATTN  ERNEST  BAUER 
OICY  ATTN  HANS  WOLFARD 
OICY  ATTN  JOEL  BENGSTON 

INTL  TEL  &  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

OICY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

11011  TORREYANA  ROAD 
P.O.  BOX  85154 
SAN  DIEGO,  CA  92138 

OICY  ATTN  J.L.  SPERLING 


JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

OICY  ATTN  DOCUMENT  LIBRARIAN 
OICY  ATTN  THOMAS  POTEMRA 
OICY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
OICY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
OICY  ATTN  DASIAC 
OICY  ATTN  WARREN  S.  KNAPP 
OICY  ATTN  WILLIAM  MCNAMARA 
OICY  ATTN  B.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 
OICY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
OICY  ATTN  DEPT  60-12 
OICY  ATTN  D.R,  CHURCHILL 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC. 

3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

OICY  ATTN  MARTIN  WALT  DEPT  52-12 
OICY  ATTN  W.L.  IMHOF  DEPT  52-12 
OICY  ATTN  RICHARD  G.  JOHNSON  DEPT  52-12 
OICY  ATTN  J.B.  CLADIS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 
OICY  ATTN  R.  HEFFNER 


M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 

OICY  ATTN  DAVID  M.  TOWLE 
OICY  ATTN  L.  LOUGHLIN 
OICY  ATTN  D.  CLARK 
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MCDONNEL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
OICY  ATTN  N.  HARRIS 
OICY  ATTN  J.  MOULE 
01CY  ATTN  GEORGE  MROZ 
01CY  ATTN  W.  OLSON 
01CY  ATTN  R.W.  HALPRIN 
01CY  ATTN  TECHNICAL  LIBRARY  SERVICES 

MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
01CY  ATTN  P.  FISCHER 
01CY  ATTN  W.F.  CREVIER 
01CY  ATTN  STEVEN  L.  GUTSCHE 
01CY  ATTN  R.  BOGUSCH 
01CY  ATTN  R.  HENDRICK 
01CY  ATTN  RALPH  KILB 
01CY  ATTN  DAVE  SOWLE 
01CY  ATTN  F.  FAJEN 
01CY  ATTN  M.  SCHEIBE 
01CY  ATTN  CONRAD  L.  LONGMIRE 
01CY  ATTN  B.  WHITE 

MISSION  RESEARCH  CORP. 

1720  RANDOLPH  ROAD,  S.E. 

ALBUQUERQUE,  NEW  MEXICO  87106 
01CY  R.  STELLINGWERF 
01CY  M.  ALME 
01CY  L.  WRIGHT 

MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

01CY  ATTN  JOHN  MORGANSTERN 
01CY  ATTN  G.  HARDING 
01CY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
01CY  ATTN  W,  HALL 
01CY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

OICY  ATTN  E.C.  FIELD,  JR. 


PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

OICY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01801 

OICY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 

OICY  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

R  &  D  ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
OICY  ATTN  FORREST  GILMORE 
OICY  ATTN  WILLIAM  B.  WRIGHT,  JR. 
OICY  ATTN  ROBERT  F.  LELEVIER 
OICY  ATTN  WILLIAM  J.  KARZAS 
OICY  ATTN  H.  ORY 
OICY  ATTN  C.  MACDONALD 
OICY  ATTN  R.  TURCO 
OICY  ATTN  L.  DeRAND 
OICY  ATTN  W.  TSAI 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
OICY  ATTN  CULLEN  CRAIN 
OICY  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

OICY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  END  AVENUE 
NEW  YORK,  NY  10023 

OICY  ATTN  VINCE  TRAPANI 


SCIENCE  APPLICATIONS,  INC. 

P.O.  BOX  2351 
LA  JOLLA,  CA  92038 

OICY  ATTN  LEWIS  M.  LINSON 
01CY  ATTN  DANIEL  A.  HAMLIN 
01CY  ATTN  E.  FRIEMAN 
01CY  ATTN  E.A.  STRAKER 
01CY  ATTN  CURTIS  A.  SMITH 
01CY  ATTN  JACK  MCDOUGALL 

SCIENCE  APPLICATIONS,  INC 
1710  GOODRIDGE  DR. 

MCLEAN,  VA  22102 
ATTN:  J.  COCKAYNE 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 

01CY  ATTN  DONALD  NEILSON 
01CY  ATTN  ALAN  BURNS 
01CY  ATTN  G.  SMITH 
01CY  ATTN  R.  TSUNODA 
01CY  ATTN  DAVID  A.  JOHNSON 
01CY  ATTN  WALTER  G.  CHESNUT 
01CY  ATTN  CHARLES  L.  RINO 
01CY  ATTN  WALTER  JAYE 
01CY  ATTN  J.  VICKREY 
01CY  ATTN  RAY  L.  LEADABRAND 
01CY  ATTN  G.  CARPENTER 
01CY  ATTN  G.  PRICE 
01CY  ATTN  R.  LIVINGSTON 
01CY  ATTN  V.  GONZALES 
01CY  ATTN  D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

01CY  ATTN  W.P.  BOQUIST 

TOYON  RESEARCH  CO. 

P.O.  Box  6890 
SANTA  BARBARA,  CA  93111 
01CY  ATTN  JOHN  ISE,  JR. 

OICY  ATTN  JOEL  GARBARINO 


TRW  DEFENSE  &  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
OICY  ATTN  R.  K.  PLEBUCH 
OICY  ATTN  S.  ALTSCHULER 
OICY  ATTN  D.  DEE 
OICY  ATTN  D/  STOCKWELL 
SNTF/1575 

VISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
OICY  ATTN  W.  REIDY 
OICY  ATTN  J.  CARPENTER 
OICY  ATTN  C.  HUMPHREY 


26 


